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Abstract 

We calculate the cross section of J/ij: inclusive production in neutrino- nucleon 
deep-inelastic scattering via the weak neutral current within the factorization for- 
malism of nonrelativistic quantum chromodynamics. Besides J/^ single production 
via the Z-gluon fusion mechanism, we also consider J plus hadron-jet associ- 
ated production. We take into account both direct production and feed-down from 
directly-produced heavier charmonia. We present theoretical predictions for the 
J/^ transverse-momentum and rapidity distributions, which can be measured in 
the CHORUS and NOMAD experiments at CERN, including conservative error es- 
timates. In order to interpret a recent CHORUS measurement of the total cross 
section, we also estimate the contribution due to J/V' prompt production via diffrac- 
tive processes using the vector-meson dominance model. 

PACS numbers: 12.38.-t, 13.60.Hb, 13.60.Le, 14.40.Gx 



1 Introduction 



Since its discovery in 1974, the J /ip meson has provided a useful laboratory for quantita- 
tive tests of quantum chromodynamics (QCD) and, in particular, of the interplay of per- 
turbative and nonperturbative phenomena. The factorization formalism of nonrelativistic 
QCD (NRQCD) ^ provides a rigorous theoretical framework for the description of heavy- 
quarkonium production and decay. This formalism implies a separation of short- distance 
coefficients, which can be calculated perturbatively as expansions in the strong-coupling 
constant Og, from long-distance matrix elements (MEs), which must be extracted from 
experiment. The relative importance of the latter can be estimated by means of velocity 
scaling rules, i.e. the MEs are predicted to scale with a definite power of the heavy-quark 
[Q) velocity v in the limit v <^ 1. In this way, the theoretical predictions are organized 
as double expansions in as and v. A crucial feature of this formalism is that it takes into 
account the complete structure of the QQ Fock space, which is spanned by the states 
n = -ft^ith definite spin S, orbital angular momentum L, total angular momentum 

J, and colour multiplicity c = 1, 8. In particular, this formalism predicts the existence of 
colour-octet (CO) processes in nature. This means that QQ pairs are produced at short 
distances in CO states and subsequently evolve into physical, colour-singlet (CS) quarko- 
nia by the nonperturbative emission of soft gluons. In the limit f — 0, the traditional 
CS model (CSM) 0| is recovered. The greatest triumph of this formalism was that it was 
able to correctly describe the cross section of inclusive charmonium hadropro duct ion 
measured in pp collisions at the Fermilab Tevatron |Q , which had turned out to be more 
than one order of magnitude in excess of the theoretical prediction based on the CSM. 

In order to convincingly establish the phenomenological significance of the CO pro- 
cesses, it is indispensable to identify them in other kinds of high-energy experiments as 
well. Studies of charmonium production in ep photoproduction, ep and uN deep-inelastic 
scattering (DIS), e~^e~ annihilation, 77 collisions, and 6-hadron decays may be found in 
the literature; see Ref. |^ and references cited therein. Here, denotes a nucleon. Fur- 
thermore, the polarization of charmonium, which also provides a sensitive probe of CO 
processes, was investigated 0,0- Until very recently, none of these studies was able to 
prove or disprove the NRQCD factorization hypothesis However, preliminary data of 
77 — > J/ip + X taken by the DELPHI Collaboration [^j at LEP2 provide first independent 
evidence for it 

In this paper, we revisit J/ip inclusive production in uN DIS. In particular, we consider 
the process u + N—^u + J/iIj + X, which is mediated via the weak neutral current (NC). 
First experimental evidence for this process was delivered two decades ago by the CERN- 
Dortmund-Heidelberg-Saclay (CDHS) Collaboration, who exposed an iron target in the 
wide-band neutrino beam produced by protons from the CERN Super Proton Synchrotron 
(SPS) and identified the J/ip mesons through their decays to pairs |jlO[- Similar 

experiments are being performed by the NuTeV Collaboration at Femilab (Experiment 
E815) |1T1 and by the CHORUS ^ and NOMAD fTSf Collaborations at CERN, using 
iron, lead, and iron targets, respectively. While NuTeV does not yet see evidence for NC 
production of J/ip mesons, CHORUS nicely confirms the CDHS measurement. On the 
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other hand, NOMAD has not yet released any results on this. 

The process v+N — > v+J/ip+X can either proceed diffractively or nondiffractively. In 
diffractive processes, the target nucleons interact with the rest of the process via exchanges 
of colourless objects with small spacelike virtualities t, and they typically stay intact 
(elastic scattering) or turn into a hadronic system of small invariant mass (diffractive 
dissociation). On the other hand, in nondiffr active processes, they participate in the 
hard scattering via their quark and gluon content and are thus destroyed (deep-inelastic 
scattering). In the case of elastic scattering off nucleons bound in compound nuclei, 
the recoiling nuclei can emerge unharmed (coherent scattering) or break up (quasielastic 
scattering) . 

The diffractive process v + N v + J /ip + N', where A^' denotes a system with little 
or no excitation, was studied in the framework of the vector-meson dominance model 
(VDM) |]I4,15,16|. The nondiffractive process u + N u + J/ip + X was treated in 



the parton model of QCD endowed with two alternative approaches: on the basis of the 
so-called Z-gluon fusion mechanism [|15],|1^,|1^ in the colour evaporation model based on 
local parton-hadron duality |]TB|; and in the GSM |T^. For reviews, see Refs. 



In all these processes, the J/ ip mesons can be produced either directly or via radiative 
or hadronic decays of heavier charmonia, such as Xcj and ip' mesons, where J = 0, 1, 2. The 
respective decay branching fractions are -B(xco J 1'^ + 7) = (0.66 ± 0.18)%, B{xci 
J/^ + 7) = (27.3 ± 1.6)%, B{xc2 -> W + 7) = (13.5 ± 1.1)%, and B{i)' J/^) + 
X) = (55 ± 5)% It has become customary to collectively denote the J/ip mesons 

produced directly or via the feed-down from heavier charmonia as prompt. If the vN 
centre-of-mass (CM) energy is sufficiently large, \fS > 2Mb ~ 10.6 GeV, then J/ip 
mesons can also originate from weak decays of b hadrons produced through the process 
u + N^u + B + B + X in diffractive or nondiffractive DIS. The corresponding decay 
branching fraction is B{B J/ip + X) = (1.16±0.10)% In fixed-target experiments, 
the requisite incident-neutrino energy E is E > (4M^ — m^) / (2m) ^ 59 GeV, where m is 
the nucleon mass. In the case of the GERN wide-band neutrino beam, only an insignificant 
fraction, namely 7%, of the total fiux satisfies this criterion [1^. Therefore and because 
of the smallness of B{B —>■ J/ip + X), J/tp production through 6-hadron decay should 
be greatly suppressed in the GHORUS and NOMAD experiments and is not considered 
here. 

In this paper, we present complete and up-to-date predictions for prompt J/ip pro- 
duction in uN NG DIS, including both the diffractive and nondiffractive contributions, 
which can be readily compared with available and expected GHORUS and NOMAD data. 
Specifically, we calculate the nondiffractive cross section of u+N — * u+J/ip+X, where X 
denotes the nucleon remnant and possibly one additional hadron jet, to lowest order (LO) 
in the framework of the NRQGD factorization formalism , and we reanalyze the diffrac- 
tive cross section of p + N — »• u + J/ip + N' in the VDM adopting the results of Refs. [0,0 
complemented with new input information. As for the nondiffractive direct contribution, 
our analysis is, in a way, complementary to the one of Ref. P2[, where e+p — > e + J/ip + X 



via photons with large virtualities was investigated: while the photon-quark coupling 
is purely vectorial, Z bosons couple to quarks mainly through axial-vector couplings. 
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the vector- coupling strengths being numerically small. However, by charge-conjugation 
invariance, the leading CSM prediction only involves vector couplings and may thus be 



obtained from Ref. [22| by appropriate substitutions [IS]. As mentioned above, the non- 
diffractive prompt CS contribution was already studied in Ref. [|l9i, two decades ago. 



However, since the authors of that paper did not present analytic results and used input 
parameters that are now obsolete, it is indispensable to repeat this analysis. 



The leading MEs of the S'-wave charmonia ^ = J/^,^' are lO"^ P^f^^l), (O^ Vs^ 



are 



3^(8) 
QXcj 



3p 



and (O^ 
and 



3p 



(8) 
J 



, with J = 0, 1, 2, while those of the P-wave charmonia x 



cj 



(1) 
J 



QXcJ 
QXcJ 



'SI 



(8) 



3p{8) 
3p{l) 
"3^(8) 



They satisfy the multiplicity relations 



3p(8) 
-I r 



(2J + 1)(C^ 
(2J+1)(C»^=« 
(2 J + 1) 



3p(l) 
3^(8)" 



which follow to LO in v from heavy-quark spin symmetry. 

At 0{a'^as), where a is Sommerfeld's fine-structure constant, the Z-gluon fusion 
mechanism [|15],|16|,|17| is realized by the CO partonic subprocesses v + g v + cc[n], 

1q(8) 3c<{8) 



where n 



3p 



(8) 
J 
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In order to enable the production of the CS Fock states 
n = 'Si^\ 'Pj^\ we need to allow for one additional gluon in the final state. The resulting 
LO CS partonic subprocesses, u + g u + cc[n] + g, are of 0{a'^al) (see Fig. |1]) [ p!9|] . 

r3p(i)i 



The cross sections of u + g ^ u + cc 



+ g suffer from infrared (IR) singularities in 



the limit of the outgoing gluon being soft. This reflects a well-known conceptual deficiency 
of the CSM 0. In the NRQCD framework these IR singularities are factorized at some 



cut-off scale and absorbed into redefinitions of the CO MEs 
the cross section of u + g ^ u + cc 'S[^^ 



QXc 



^ ) that multiply 



In turn, these CO MEs become scale dependent 
at 0{as) and satisfy appropriate evolution equations. The factorization is conveniently 
performed in dimensional regularization p4|. In Ref. |T^, these IR singularities were not 
encountered because the hadronic recoil energy Eh was required to exceed some minimum 
value -E™™. However, a rapid rise of cross section with decreasing value of i?™™ was 
observed. 

In the same order, we also have the CO partonic subprocesses u + a v + cc[n] + a, 
where a = q,q, g, with q = u,d, s, and n = ^Sq^\ 'S[^\ 'Pf'^ (see Fig. p. The cross sections 
of the latter are plagued by collinear singularities in the limit of vanishing cc transverse 
momentum in the Z*N CM frame. Here, Z* denotes the virtual Z boson. These 
singularities could be avoided by introducing an appropriate cut on p^. In a full next- 
to- leading-order (NLO) analysis, they would be factorized at mass scale /i/ and absorbed 
into the bare parton density functions (PDFs) of the nucleon so as to renormalize the 
latter. 

In this paper, we provide analytical results for the cross sections of all the 2^2 and 
2^3 partonic subprocesses enumerated above to LO within the NRQCD factorization 
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formalism In the 2 — >■ 2 case, we found agreement with Eqs. (8) and (9) of Ref. [p3[]f], 
so that there is no need to hst our formulas for the squared transition-matrix elements. 
We then calculate the inclusive cross section of prompt J / ip production in nondiffr active 
pN NC DIS under CHORUS and NOMAD kinematic conditions to LO in NRQCD and 
the CSM. This involves the 2^2 CO processes of J/V', Xcj, and ip' production and 
the 2 ^ 3 CS processes of J/ip and ip' production. As explained above, the 2-^3 
CS processes of Xcj production lie beyond the scope of the CSM. We also consider the 
distributions in the J /tl) transverse momentum px and rapidity y in the laboratory frame, 
which should be experimentally accessible. 



Table 1: Relative importance, measured in powers of and v"^ ^ as, of the various 
nondiffr active J/ip production channels according to the NRQCD counting rules. The 
entries are normalized to the direct 2^3 CS contribution, which is of O (a'^a'^v^). 





J/i^, i^' Xcj 


2 ^ 2 CO 
2 ^ 3 CS 
2^ 3 CO 


v'^/as^as v^/as^l 
1 v"^ as 
~ fts f ^ ~ as 



The relative importance, measured in powers of as and ~ a^, of the various non- 
diffractive J / ip production channels according to the NRQCD counting rules is represented 
in Table 0. The entries are normalized to the direct 2 — »• 3 CS contribution, which is of 
O {a'^a^v^). In addition, the feed-down channels are suppressed by the respective branch- 
ing fractions. In fixed-target experiments, which are typically performed at low values 
of [10,11,12,13|, phase-space suppression comes in as another limiting factor. This 



affects the 2-^3 processes more severely than the 2 — >■ 2 ones. In fact, as will become 
apparent in Section |^, this is the reason why, in the case of direct J/ ip and ip' production, 
the 2 — >■ 2 CO contributions greatly exceed the 2 — >■ 3 CS ones, although they are formally 
suppressed according to Table |I|. If events of nondiffr active direct J/ip and ip' produc- 
tion could be identified experimentally, this would lend itself to a powerful discriminator 
between the CSM and NRQCD. Obviously, the 2^3 CO contributions, which are not 
included in our numerical analysis for the reasons explained above, would only provide 
minor corrections to their 2 — 2 counterparts, which make up the bulk of the NRQCD 
prediction. 

This paper is organized as follows. In Section we present, in analytic form, the cross 
sections of the partonic subprocesses enumerated above and explain how to calculate from 
them the total cross section of prompt J/ip production in nondiffractive uN NC DIS as 
well as its pt and y distributions. Lengthy expressions are relegated to the Appendix. 

^There is a typographical error in the last equation of Eq. (9): (Q^ + 4m^)'^ in the denominator should 
be replaced by (Q^ + 4rn^)^. 
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For the reader's convenience, we also recall the cross-section formulas for prompt J/ip 
production in diffractive vN NC DIS. In Section we present our numerical results and 
compare them with recent CHORUS data. Our conclusions are summarized in Section |^. 



2 Analytic results 

In this section, we present our analytic results for the cross sections of v+N u+H+j+X 
in nondiffr active DIS, which proceeds through the 2 — 3 partonic subprocesses u + a ^ 
V + cc[n] + a mentioned in Section |l|. Here, = J/'?/', Xcj, V'') the hadron jet j arises from 
the fragmentation of the additional final-state parton a, and X is the nucleon remnant. 
In our numerical analysis, we only consider the CS processes of Jj^l) and ■?/'' production 
and integrate their cross sections over all kinematically allowed values of the j three- 
momentum. The CS processes of Xcj production and the CO processes are provided for 
future applications. At the end of this section, we also present formulas for the cross 
sections of the processes v^X-^v^H^X and v^X—^v^H^ X' . 

We work at LO in the parton model of QCD with = 3 active quark flavours and 
employ the NRQCD factorization formalism to describe the formation of the H meson. 
We start by defining the kinematics. As indicated in Fig. ^, we denote the four-momenta 
of the incoming neutrino and nucleon and the outgoing neutrino, H meson, and hadron 
jet by fc, P, fc', 'Pill and p', respectively. The parton struck by the virtual Z boson carries 
four-momentum p = xP . We neglect the masses of the nucleon and the light quarks, call 
the one of the H meson M, and take the charm-quark mass to be mc = M/2. In our 
approximation, the nucleon remnant X has zero invariant mass, M\ = [P — •pf' — 0- 
The CM energy square of the vX collision is 5 = (A; + P)^. The virtual Z boson has 
four-momentum q = k — k' . As usual, we define = — > 0, y = q ■ P/k ■ P, and 
the inelasticity variable z = pn ■ P/q ■ P. In the nucleon rest frame, y and z measure the 
relative neutrino energy loss and the fraction of the Z* energy transferred to the H meson, 
respectively. The Z*X CM energy square is W"^ = {q + P)"^ = yS — Q"^. The system X' 
consisting of j and X has invariant mass square M|, = (g + P — puY = (1 ~ — z)S. 
As usual, we define the partonic Mandelstam variables as s = {q + pY = xyS — Q^, 
t = {q — PhY = —xy{l — z)S, and u = {p — pnY = M"^ — xyzS. By four-momentum 
conservation, we have s + t + M = - Q2 in the Z*X CM frame, the H meson has 
transverse momentum and rapidity 



t(su + Q2M2) 
1 s{M^-u) IW"^ 



Pt = ' . , o2 (2) 



y* =-\n—. ^ — ^ h-ln , (3) 

2 s (m2 - i) + Q2M2 2s' ^ ^ 

respectively. Here and in the following, we denote the quantities referring to the Z*X 
CM frame by an asterisk. The second term on the right-hand side of Eq. (j^) originates 
from the Lorentz boost from the Z*a CM frame to the Z*X one. Here, y^ is taken to be 
positive in the direction of the three-momentum of the virtual Z boson. 
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In the parton model, the nucleoli is characterized by its PDFs fa/N{x, ^f), and, at LO, 
an outgoing parton may be identified with a hadron jet. Thus, we have 

da{u + N ^ u + H + j + X) = f dx"^ fa/N^x, ^f)da{u + a ^ u + H + a), (4) 

where a = u,u, d, d, s, s, g. Furthermore, according to the NRQCD factorization formalism 
we have 

da{iy + a ^ u + H + a) =J2 (^'^M) ^^^(z/ + a ^ u + cc[n] + a), (5) 

n 

where, to LO in n = ^S[^\^sli^\^S[^\^PP for H = J/^,^' and n = ^P^^ ^^f ^ for 
H = Xcj. 

Decomposing the transition-matrix element of the partonic subprocess u + a —>■ u + 
cc[n] + a into a leptonic part. 



where g = Mz, with Gp being Fermi's constant, and a hadronic one, T^{Z* + a — ^ 

cc[n\ +a), from which the virtual Z-boson leg is amputated, we can write its cross section 
as 

da{. + + cc[n] + a) = ^ ^ IJ^W (^^'^ ^'^'^^) 

xdVS^{k + p-k\pH,p), (7) 

where Ng = Ng = Nc = 3 and Ng = {N^ — 1) are the colour multiplicities of the 
partons a and the hadronic tensor H'^" is obtained by summing the absolute square of 
T'^{Z*+a — >■ cc[ra]+a) over the spin and colour states of the incoming and outgoing partons 
a. Here and in the following, we employ the Lorentz-invariant phase-space measure 

dPS^{p-p„ . . . = (2vr)^5W (p - ±p) n "^'^^ 



i=l 



,0- 



The first factor on the right-hand side of Eq. (|^ stems from the fiux and the second one 
from the average over the spin and colour states of the incoming parton. Integrating over 
the azimuthal angle of the outgoing neutrino, we may simplify Eq. (|^) to become 

da{u + a^ u + cc[n] + a) = 77^ tttt ('t) L'^'H^^— ^ ^ dV^2{q + P]PH,p'), 



2xS 2Na UnJ y (Q2 + jvf|)' 

(9) 

where 

^ i + i^~y)\ fu. _ ^i-y) ^fu. ^ ^2 _ yyj^^ (lo) 



with 



tj. — 




^^^^ _ 




— 








-V 



is the leptonic tensor. Here, we adopt the convention 60123 = 1- The cross section of 
17 + a — i> z7 + cc[n] + a emerges from Eq. (^) through crossing symmetry, by flipping the 
sign of the last term on the right-hand side of Eq. (]T0|). In the following, the symbol v 
collectively denotes neutrinos and antineutrinos. 

We evaluate the cross sections of the relevant partonic subprocesses p+a — > i'+cc[n]+a 
from Eq. (^) applying the covariant-projector method of Ref. ||25|. Our results can be 
written in the form 



a 



dy dQ"^ di 



[u + a ^ u + cc[n] + a) = Fa[n] 



Ta\n] 



4(1-?/) 



Lain] 



+ {2-y)A4n] 



(12) 



where -FaM, Ta[n\, La[n\, and Aa[n\ are functions of s, t, -u, and Q^, which are listed 
in the Appendix. We combined the results proportional to the CO MEs (O"^ 



and (O^-J 



^S[^^ ^ exploiting the multiplicity relations of Eq. (|TD. The Ta[n] and La[ 
functions involve terms proportional to ViVj or aiaj, while the Aa[n\ functions involve 
terms proportional to Vittj, where i,j = q,c. Here, Vq = — 2eg sin^ 9w and 
the Zqq vector and axial-vector couplings, respectively, where is the third component 
of weak isospin of the left-handed component of quark g, is the fractional electric charge 
of the latter, and 6^^ is the weak mixing angle. We recover our result for the cross section 
ofe + a^e + J/ip + a, given in Eq. (13) of Ref. P, by substituting in Eq. (0) g = e, 
Vq = Cq, Vc = Cc, = Oc = 0, and Mz = 0, where e = A/ivra is the electron-charge 
magnitude. 

For the CO Fock states n = ^S^^\^S[^\^pP , the cross sections of Eq. (|12D exhibit 
coUinear singularities in the limit t —>■ 0. According to the factorization theorem, the 
limiting expressions must coincide with the respective u + g —>■ v + cc[n\ cross sections [p^ 
multiplied by the spacelike a ^ g splitting functions. This provides another nontrivial 
check for our results. 



3p(8) 



n\ 



It is interesting to observe that the cross sections oi v + g ^ v + cc 



in the limit ?/ — and ■ 
Inserting Eq. (|^) in Eq. 

over X and t, we obtain 



3p 



(1) 
J 



+ g vanish 



0. 



and including the maximum boundaries of the integrations 



dy dQ"^ 



[u + N^u + H + j+X) 



dx 



(Q2+M2)/(yS) 



(s+Q2)(s_M2)/s 



dt 
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a 



dy dQ"^ di 



V + a — z/ + cc[n] + a), 



(13) 



where {d? a / dy dQ'^ dt^ {u + a ^ u + cc[n] + a) is given by Eq. (|T2D. The kinematically 
allowed ranges of 5", y, and are S > M^, M'^/S < y < 1, and < Q"^ < yS - M^, 
respectively. We then evaluate the cross section of prompt J/ip production as 

da{u + N + ( J/V^)prompt +j + X) = da{v + N^v + J/i) + i + X) 



+ E B{xcj ^ W + l)d(y{y + N ^v + Xcj + J+X) 



j=o 



+ B{^' Jj^ + X)do[v + ^ z/ + V^' + J + X). 



(14) 



The distributions in y and can be evaluated from Eq. (13) as it stands. It is 
also straightforward to obtain the distributions in and y*^^ given in Eqs. (H) and (H), 
respectively, by accordingly redefining and reordering the integration variables in Eq. (pTS]). 
The distribution in the J/^ azimuthal angle 0* in the Z'^X CM frame is constant. 

The evaluation of the distributions in the J/V' transverse momentum p^, rapidity yu-, 
and azimuthal angle in the laboratory frame is somewhat more involved. Choosing a 
suitable coordinate system in the Z*X CM frame, we have 



{k 



{PhY 



S-Q 



2W 



2W 



sin 


Xcosip* J 






2W 



2W 



/ cosh y"^ \ 
cos (j)* 
sin 0* 
V my sinh y^ j 



i 1 \ 
sin^ 



Vcos^V 




V-1/ 



(15) 



whereQ 



COS'?/''* 



2SW'^ 



(5-Q2)(H^2 + g2) 



1, 



cos^ 



1 



(5- iy2)(|y2^_Q2) 



m;^ = V + br) • (16) 
On the other hand, in the laboratory frame, which coincides with the nucleon rest frame. 



^The formula for cos d* given in Ref. p2[ contains a sign error, which was absent in the preprint version 
of that paper. 
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we have 



/1\ 






E' 



\ 



Ph 



where 

E 

cos ip 



2m ' 



-qsinip 


V qcosip ) 

I rriT cosh yn \ 
Pt cos 
Pt sin 

V mT sinh i/h J 



S-W^-Q^ 



I 1 
sin Q 


Vcos^y 






(17) 



2m 



cos 6* 



1 - 



1 
1 



2m ' 
2m 



+ 



Q'^m 
S — m^ 



mx = JM^ + p^. 



(18) 



Notice that yn is taken to be positive in the direction of the three-momentum of the 
incoming neutrino. Without loss of generahty, we may require that < ip^iip < it, 
for, otherwise, we can achieve this by rotating the respective coordinate systems by 180° 
around the z axis. We can then evaluate pt, yn, and from p^, y'^, and 0* as 



Pt 
Vh 



(p^) + A{A-2p^cos(f)*), 
, , iW^ + Q^)m*T 1, S 



cos I 



P^ COS ( 



-A 



Pt 



where 



A 



-mTexpiy^) 



(19) 
(20) 

(21) 
(22) 



1 + cos Ip* V 

The third term on the right-hand side of Eq. (pOD stems from the Lorentz boost from the 
uN CM frame to the laboratory one. Since px, yn, and depend on 0*, the integra- 
tion over 0* in Eq. (|l^) is no longer trivial, and we need to insert the symbolic factor 
(l/27r) Jq"^ d(p* on the right-hand side of that equation. 
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For future applications, we also present compact formulas that allow us to determine 
p^, y^, and 0*, once px, Vh, and are given. In fact, Eqs. ([T9|) and (|2l|) can be straight- 
forwardly inverted by observing that the quantity A defined in Eq. (|22|) can be expressed 



in terms of m-r and yn by substituting Eq. (poD , the result being 



Having obtained p^, we can then evaluate from Eq. (IzO]). For the reader's convenience, 
we collect the relevant formulas here: 

Pt = \Jpt + A{^ + '^Pt COS0), 

^H = l/H + ln ^^,^^,^^^ + -ln-, 

cos0- = ^"^"\^ + ^ . (24) 
Pt 

If the four-momentum p/^ + p' + P — p of the hadronic final state H + j + X can 
be measured in the laboratory frame, then the Z* four-momentum can be evaluated as 
Q = [Ph + p' + P — p) — P- Since the direction of the incident neutrino beam is given 
by the experimental set-up, the four-momenta k and k' of the incident and scattered 
neutrinos can thus be reconstructed using Eqs. {^n\j and (plSf). Then, the Z*N CM frame 
is well-defined, and p^, y^, and (p* can be determined. However, if the hadronic final 
state cannot be fully detected, then one can still measure pt, Vh, and in the laboratory 
frame. 

We now turn to the cross section of u + N u + H + X in nondiffractive DIS, 
which proceeds through the 2 — > 2 partonic subprocesses u + g ^ u + cc[n\ mentioned in 
Section |I]. We have 

-^^(i/ + N^iy + H + X) = — TV / dx fa/Nix, Uf) 

^Y.{o"[n])K{y.Q^). (25) 



where, to LO mv,n = ^S'i\^Sf\ ^Pf for H = J/^, i)' and n = ^ for H = Xcj- The 
functions hn{y,Q^) may be found in Eq. (9) of Ref. [^], which contains a typographical 
error, as pointed out in Section |I|. The kinematically allowed ranges of S and are 
S > M'^ and < < S - M^, respectively. The kinematics oiu + N^u + H + X 
emerges from that of u + N ^ u + H + j + Xhj nullifying p'. This leads to the simplified 
relations xy = {Q^ + M^)/S, z = 1, Mx' = 0, s = M^, t = 0, m = -Q^, p1^ = 0, 



y^ = ln(iy/M), Pt = y {I — y)Q^, yn = \n-{yS/mmT), and cos0 = —1. As before, we 
have W"^ = yS — Q"^ . In turn, and x may be expressed in terms of pr and yn, as 

g2 ^ Pt ^ ^ mTexp(-yff)/m - M"^ / S 

1 ~ mniTQ'X-viyH)/ S' 1 — mniT e^Y>{y h) / S ' 
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which allows us to conveniently evaluate the pt and yn distributions from Eq. (p5|). 

In the remainder of this section, we recall the analysis [|l5| , p!6| of the total cross section 
ofi' + N—>-i' + H + N', with H = J/ip, Xci, 'ip', in diffractive DIS. We first consider the 
case H = J/ip. The differential cross section ofu + N—yi' + H + N' may be be obtained 
from the one of fi + N —>■ fi + H + N' in diffractive DIS via a virtual photon 7* by adjusting 
the electroweak couplings and the propagator, as 



_fa_ 

dydQ"^ 



(u + N ^ u + H + N') 



(Q2 + M|) 



dy dQ' 



ifi + N ^ fi + H + N'). (27) 



The latter can be decomposed as |26 



dy dQ' 



a y-Q^/{2Em) 
2^ Q2{y^ + QyE^: 



+ 



2(1 -y) 



2^2 



l + {l-yf + 



2^2 



:28) 



where (Jt{i^, Q^) and o"l(z/, Q'^) are the transverse and longitudinal parts of the diffractive 
cross section of '-/* + N ^ H + N' in virtual photoproduction, respectively. Here, E = 
{S — m2)/(2m) and u = q ■ P/m = yE = g° are the incoming-lepton and virtual-photon 
energies in the nucleon rest frame, respectively, is the lepton-nucleon CM energy, and 
the muon mass is neglected. At small values of Q^, the contribution from longitudinal 
photons is suppressed, and we may put aL^UjQ'^) = 0. On the other hand, according 
to the VDM, the dependence of axiy^Q"^) is approximately described by the familiar 
dipole suppression factor, viz. 



(1 + Q2/M2)2^ 



(29) 



where ct(z/) is the diffractive cross section of 7 + — > if + A^' in real photoproduction and 
My is a phenomenological mass parameter. Fits to experimental data of 7* + A^ H + N' 
suggest that My = (3.2 ± 0.6) GeV ||2^. The measured z/ dependence of (t(z/) is well 



described by the phenomenological parameterization |T3|,|TB|,| 



0- z/ 



Aexp 



-B 



(30) 



with A = 20 nb, 5 = 45 GeV, and C = 6 GeV [0, which is valid for v>C. 

We now generalize our considerations to also include the cases H = Xci^i^'- To this 
end, we recall that the strengths of the Z-boson vector and axial-vector couplings to the 



S- and P-wave charmonia are proportional to Vcfv and flc/yi; where [16 



^2 _ 3 \Rs 



mi 



9 \R'n 



(31) 
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respectively. Here, Rs{r) and Rp{r) denote the radial wave functions of the S- and 
P-wave charmonia, respectively. Furthermore, we have [|l| 



\R, 



\R!Am 



o 



2tx 
"9 



3p(l) 



2J+ 1 



(32) 



Assuming that a{il)' + N i)' + N') ^ a{xci + N Xci + N') ^ a{J/i) + N J/^ + N'), 
we thus conclude that dV(z/ + ^ z/ + + N')/{dydQ^) and rfV(i/ + ^ z/ + Xci + 
N')/{dydQ'^) may be obtained from rfV(z/ + A^ z/ + J/?/' + N')/{dydQ^) through 
multiplication with the constant factors 











-35(1) 


1) 





3p(l)] 


> 




"35(1)" 


) 



Rm 



R = ^ A 



respectively. Thus, the effective direct-to-prompt conversion factor reads 

R{ji^)„. = 1 + Rx.^Bixci ^ + 7) + i?v"5(^' ^ + ^). 

The kinematically allowed ranges of E, Q^, and y are 

M 



(33) 



(34) 



E>M{1 + 



2m 



o<g2 < 



2E[2Em - M(M + 2m)] 



2P + m 



Q2 + M(M + 2m) 
2Em 



<y<l 



4E^'' 



respectively. Here, the A^' mass is taken to be m. In addition, Eq. 
y > C/E. 



(35) 

implies the cut 



3 Numerical results 

We are now in a position to present our numerical results. We first describe our theoretical 
input and the kinematic conditions. We use m = 0.938 GeV, rric = (1.5 ± 0.1) GeV 
III, Mv = (3.2 ± 0.6) GeV Mw = 80.419 GeV, Mz = 91.1882 GeV, sin^^^ = 
1 - Mlr/Ml = 0.22225, Gp = 1.16639 x 10"^ GeV'^, a = 1/137.036, and the LO 
formula for a'r^\jji) with nj = 3 active quark flavours [E^. The CHORUS Collaboration 
1^ chose the target material to be lead (Pb), which in average consists of A = 207.2 
nucleons, Z = 82 of them being protons. The NOMAD Collaboration uses iron 
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(Fe), with A = 55.854 and Z = 26. Leaving aside nuclear corrections and appealing to 
strong-isospin symmetry, the effective nucleon PDFs may be approximated by 



1 r 



fa/N{x, /i/) = ^ Zfa/p{x, /i/) + {A- Z)fb/p{x, /i/) 



fc/N{x,fif) = fc/p{x,fif), 



(36) 



where (a, b) = {u, d), (u, d), {d, u), {d, u) and c = s, s, g. We present the cross sections per 
nucleon, rather than per nucleus. The dependence of the cross section per nucleus, cta, on 
the number of nucleons. A, was studied experimentally for J / ip inclusive photoproduction 

^H]. The New Muon Collaboration (NMC) H, who per- 



31| , |32| and muopro duct ion 



formed the most recent of these experiments, was able to reconcile the seemingly inconsis- 
tent results of the previous experiments ||3ll , |32| , |33| by imposing appropriate cuts on z and 
Pt- Separating signal events of elastic (coherent) and deep-inelastic (incoherent) scatter- 
ing, the A dependences of the respective cross sections cr^ were found to differ significantly 
||3^ , |3^ . Assuming that a a is related to the cross section per free nucleon (hydrogen), ai, 
by the simple power law a a = A°'ai [32,^], with a being an A-independent number, we 
extract the values ttcoh = 1.19±0.02 and a^^ = 1.05±0.03 from the NMC results for coher- 
ent and incoherent scattering, respectively. For lead (iron) targets, this implies that the 
elastic and deep-inelastic cross sections per (bound) nucleon, a a/ A, are enhanced relative 
to the respective values of ai by factors of 2.75+J;g (2.15+J;;5) and 1.31^S:2o 
respectively. We note in passing that the above value of a^^ is in same ball park as those 
recently found for inclusive large-p^ prompt-photon and neutral-pion hadroproduction 
in fixed-target scattering, namely 1.04 and 1.08, respectively ||35|. We assume that the 
power laws of J/ip muoproduction, which dominantly proceeds via the electromagnetic 
current, carry over to uN NC DIS and employ them to determine the nuclear corrections. 
As for the proton PDFs, we employ the LO set by Martin, Roberts, Stirling, and Thorne 
(MRST98LO) [Q, with asymptotic scale parameter A^^^ = 174 MeV, as our default and 
the LO set by the CTEQ Collaboration (CTEQ5L) |3, with A^^) = 192 MeV, for com- 
parison. The corresponding values of A*^^-* are 204 MeV and 224 MeV, respectively. We 
choose the renormalization and factorization scales to be /ij = ^i\/Q^ + M'^, with i = r, f, 
respectively, and independently vary the scale parameters and ^/ between 1/2 and 2 
about the default value 1. We adopt the NRQCD MEs from Table I of Ref. [0, which 
were determined there from the measured leptonic annihilation rates of the J/ip and ip' 
mesons and fits to the Tevatron data 0]. As for the J/ip and ^p' mesons, the fit results 



for (O^ 



lc.(8) 
On 



and ( 



(8) 



are strongly correlated, so that the linear combinations 



Mf = (O 



lc.(8) 
On 



mi 



3p{8) 



(37) 



with suitable values of r, are quoted. Since Eq. ( |T3|) is sensitive to different linear combi- 



nations of (O"^ 
n^pMf and 



15<(8) 



3n.(8) 



than appear in Eq. (|37|) , we write (O'^ 



1q(8) 



and (^O^ 

= (1 — K^) (ml/r) Mf and independently vary kjj^ and k^i be- 
tween and 1 around the default value 1/2. The wide-band neutrino beam utilized in the 
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CHORUS and NOMAD experiments mainly consists of neutrinos and contains I/^, z/g, 
and Vf. admixtures of about 6%, 1%, and below 1%, respectively. The average z/^ energy is 
27 GeV. Using accurate parameterizations for the energy spectra (f)^{E^) of the individual 
neutrino species z/ |T^, we calculate the spectrum-averaged cross section as 

where u = z/g, z7e, t'/i, J^^t- In order to conform with the CHORUS analysis, we apply the 
acceptance cut 20 < E^, < 200 GeV. 

In order to estimate the theoretical uncertainties in our predictions, we vary the un- 
physical parameters ^r, ^f, f^j/ip, and k^i as indicated above, take into account the exper- 
imental errors on Ocoh, c^m, f^c, My, the relevant charmonium decay branching fractions, 
and the default MEs, and switch from our default PDF set to the CTEQ5L one, prop- 
erly adjusting A'^'^^ and the MEs. We then combine the individual shifts in quadrature, 
allowing for the upper and lower half-errors to be different. 

In Figs. II and ||, we respectively present the pt and yj/^ distributions of prompt J/ip 
production in nondiffractive uN NC DIS under CHORUS experimental conditions pre- 
dicted to LO in NRQCD. Specifically, the 2-^2 CO contributions (upper histograms) are 
compared with the 2-^3 CS ones (lower histograms). In each case, the theoretical errors, 
evaluated as explained above, are indicated by the hatched areas. Since the contributing 
partonic cross sections are all gluon initiated, our our central predictions can be rendered 
appropriate for NOMAD by simply adjusting the nuclear correction factor, by multipli- 
cation with 0.94. We observe that the shapes of the CO and CS distributions are very 
similar. However, the CO distributions greatly exceed the CS ones in normalization. In 
fact, the integrated CO and CS cross sections are 3.0 x 10~^ fb and 1.0 I^q^ x 10""^ fb, 
respectively, the ratio of the central values being 30. The bulk parts, 72% in the CO case 
and 78% in the CS case, are due to direct J/ip production. In the first case, the fractions 
due to feed-down from Xcj and ip' mesons are each 14%, while, in the CS case, there is 
only feed down from ip' mesons. Therefore, at first sight, the suppression of the CSM 
prediction seems to contradict the naive expectation based on the NRQCD counting rules 
(see Table |I]). However, detailed inspection reveals that this suppression can be attributed 
to the smallness of the effective value of v^, which affects the 2 — > 3 phase space of the 
CS process more severely than the 2^2 one of the CO processes. For E = 27 GeV, we 
just have ^/S = 7.2 GeV. For a monoenergetic neutrino beam with E = A8 TeV, so that 
VS = 300 GeV as for ep collisions at DESY HERA, the CO to CS ratio becomes 4.5. 
This should be compared with the corresponding analysis for direct J/ijj production in 
nondiffractive ep DIS at HERA, with > A GeV^ which yields a CO to CS ratio of 7.8 
2|. 

We find the total cross section of prompt J/ip production in diffractive uN NC DIS 
to be 3.3 1*1 J^'^ X 10~^ fb. The theoretical uncertainty quoted here is merely parametric; it 
neither includes the experimental error on cr{i>) as parameterized by Eq. (^Dj), nor does it 
reflect the unavoidable lack of rigour of the underlying model assumptions. The fractions 
due to feed-down from Xci and ip' mesons amount to 19% and 17%, respectively. The 
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corresponding prediction for NOMAD is obtained through multiphcation with the scahng 
factor 0.80, which accounts for the difference in the nuclear correction factor. 

The CHORUS Collaboration measured the spectrum-averaged total cross section per 
nucleon of z/ + ^ z/ + J/^ + X in DIS to be (6.3 ± 3.0) x 10"^ fb. In Fig. |, this 
experimental data point is shown together with our theoretical predictions for the CO, CS, 
and diffractive contributions. It should be compared with their sum, 6.4 x 10~^ fb. As for 
the central values, the measurement exceeds the combined prediction by a factor of 10. We 
have to bear in mind that our theoretical prediction for the nondiffr active cross section is of 
LO in as and v. The NLO corrections are expected to lead to a substantial enhancement. 
It would be desirable if the CHORUS Collaboration were able to substantially increase 
their data sample and to exclude the diffractive regime by imposing suitable acceptance 
cuts on 2, p^, or M^', and if the NOMAD and NuTeV Collaborations could come up with 
similar measurements. 



4 Conclusions 

We provided, in analytic form, the cross sections of the partonic subprocesses v + a ^ 
V + cc{n\ + a, where a = g, g, ^ and n = '^S^l\ ^pf\ ^S^^\ ^S?\ ^Pf\ to LO in the NRQCD 
factorization formalism [|l|. We also confirmed previous results for the cross sections of 
the partonic subprocesses i' + g v + cc[n], where n = ^Sq^\'^Si^\'^Pj'^ [^. Using 
these results, we then studied the cross section of prompt J/ip inclusive production in 
nondiffr active uN NC DIS to LO in NRQCD and the CSM. We presented the cross section 
per nucleon bound in lead and averaged it over the effective energy spectrum of the wide- 
band neutrino beam employed in the CHORUS and NOMAD [|13] experiments. Apart 



from the total cross section, we also considered the px and yj/^ distributions. The cross 
sections of the 2 — >■ 3 CO partonic subprocesses, which did not enter our phenomenological 
study, can be used in the future, at LO, in connection with an appropriate acceptance 
cut on z, p^, or Mx' to exclude the collinear singularities in the limit t or, at 
NLO, as an essential ingredient for the real radiative corrections to the inclusive cross 
section of u + N—^u + J/iIj + X. A similar comment applies to the 2^3 CS partonic 
subprocesses of Xcj production, where the singularities are of IR type. Repeating the 
analysis of Refs. |TB|,|TB[ with up-to-date input, we also evaluated the total cross section 
of prompt J/ip inclusive production in diffractive pN NC DIS. 

As for nondiffr active direct J/ip production, from the NRQCD counting rules, one 
expects the CO contribution, which arises from 2 — 2 partonic subprocesses, to be 
slightly suppressed, by a factor of /as ~ as relative to the CS one, which is generated 
by 2 — > 3 partonic subprocesses. However, our analysis revealed that, under CHORUS 
kinematic conditions, the CS contribution is approximately 30 times smaller than the CO 
one. We demonstrated that this may be attributed to the smallness of the effective value 
of \fS, which constrains the 2 — > 3 phase space more severely than the 2^2 one. 

As far as central values are concerned, the CHORUS measurement of the spectrum- 
averaged total cross section per nucleon oiv + N^v + J/%l) + Xui DIS overshoots our 



16 



combined prediction by a factor of 10. However, this comparison does not yet permit any 
meaningful conclusions concerning the validity of the NRQCD factorization formalism 
|l|]. On the one hand, our theoretical prediction for the nondiffr active contribution should 
be substantially increased by the inclusion of NLO corrections, which are presently still 
unknown. On the other hand, the experimental error is still rather sizeable, and the 
diffractive events have not been separated from the experimental data set. 
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A Partonic cross sections 

In this Appendix, we present analytic expressions for the coefficients -FaM, ^al'^], -^aN, 
and Aa[n] appearing in Eq. (|1^). In order to compactify the expressions, it is useful to 
introduce the Lorentz invariants s = 2q-p, t = —2p-p', and u = —2q-p', which are related 
to the partonic Mandelstam variables by s = s + Q^, t = and u = u + Q"^, respectively. 



i^{y) + q{q) 



v{v) + q{q) 



Z/(z7) + CC 
Z/(z7) + CC 
Z/(z7) + CC 

F-- 



3p(l) 



1^(8) 



F = T 
+ q{q): 
F = T 



A = 0. 



L = A = 0. 



+ q{q): 



u 



T = 2QV + 2Q''st{s + u) + s* + 



L 

A 



-2QH{QH 



su) 



0. 



(A.l) 



(A.2) 



(A.3) 
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u(v) + q(q) ^ u(v) + cc ^1 + q{q): 



T = -(aj + v^^)QH{s + uy{2QH\2s + t) + 2Q^s[s^u - st{3t - 2u) - 2t^] + Qh^[s\t 

- 2u) + 2s(t2 - 4tM - u^) + t{2t^ - 2tu - u^)\ + s^u{s^ + 2st + + 2tu + 

+ 4agacQ^t(s + u){(f - s){Q^ - u){QH{s + 2t) + Q^[s^ - s^{t - 2u) - bst-^ - 2t^{t 
+ u)] - Q'^s[s^ + s^(t + 5m) - s{2t^ - btu - 3m^) - t{2t^ + tu- u"^)] + s'^u[2s'^ + s{At 
+ 3u) + + 3tu + u^]} - 2alQ\Q^ - sf{Q'^ - uf{2QH'^ + 2QH[s^ - s{3t - u) - 3t{t 
+ u)] - Q%s^ + 6sH + s^{2f + mu + v^) - 2st{Af + tu- 3u'^) - Af{t + uf] + 
+ s^(5i + u) + s\9>t^ + I2tu + u^) + s\At^ + m'^u + mu^ + u^) + 4stM(t + uf}, 
L = 2(a5 + vl)Q\Q^ - s)H\s + tf{s + uf - iaga^Q^Q^ - sf{Q^ - u)t{s + u)[Q%s 

+ 2t) - Q'^{s^ + 2sH + 3st^ + 2t^ + 2i^M) + s'' + s^(2i + u) + sH'^ - stu(t + u)] + 2a^(g2 

- s)^(Q^ - uf{2Qh'^ - 2Qh[s^ + s{3t - u) + 3t{t + u)] + Q^[s^ + AsH + s^(8t^ - 2tu 
+ u^) + 2st{Af + tu- 3u'^) + Af{t + uf] - 2Q'^s[s^ + s^{2t + ii) + s^{2t^ + v?) + s{f 

- 2t^u - 2tu^ + u^) - 2tu{t + uf] + + t + uf{s^ + u^)}, 
A = ±2vqQ'^st{s + u){aq{s + 'u)[2gSt - 2QH{2s^ + st - tu) 

+ - 2u) + 2s(i2 + 14^) - tu{2t - u)) + s^u{s -u){s + 2t + u)\ 

- 2ac(Q^ - s)(g^ - ii)[g^i - Q^is^ + 3si + i^) 

+ Q'^{s^ + s^(3i + li) + s{2t^ + tu + u^) - tu{t + u)) - su^{s + 1 + u)]}, (A.4) 
where the plus (minus) sign refers to a u + q or 17 +q + q or V + q) initial state. 



3p{8) 



+ q{q): 



37rM3(Q2 + M|)2s4t(s + 



■u 



i4' 



T = 8g^t[2s^ + st + t{2t + ?7,)] - 2g^[4s^ + l2sH + s^(19t^ + Stu + 4^^) + 2st{<6t^ 
+ tu- 2u^) + t\8t^ + 12tu + 7u^)] + 2Qh[6s^ + s^{7t + 6u) + s^{4:t^ + 3tu 
+ 6m^) - sM(8t^ + 3tM - 6m^) - tM(8t^ + I2tu + 7m^)] - s^(s + m)[7s^ + s^(12t 
+ 7m) + s(8t2 + 16tM + 7u^) + M(8t2 + I2tu + 7^^)], 

L = 8g^t[s^ - st - t(2t + u)] - 2Q^[2s^ + 4s^t + s'^{^t^ + 8tM + 277,^) - 2st(6t^ 
+ tM - 2u^) - t^(8t2 + 12tM + lu^)] + 2g^s(s + m)[2s^ + 2s^t + si^t^ + btu 
+ 2^2) + t(8t2 + i2tM + Tm^)], 

A = 0. (A.5) 



i/(i/) + 51 ^ i^(i^) + cc 



-2gW,vlQ^ 



277rM(g2 + M|)2s4(s + i)2(s + ^^)2(^ + ^^)2 ' 
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T = iQH'^is'^ + t^) - 2QH[s^{?,t - 2u) + ?,sH{t + -u) + 2st^{t - u) + 2t^{t + u)] 
+ 2Q^s[s^{t - uf - 2sHu{t + u)- st^u{2t - u) - 2t^u{t + u)] - 2s^[s^{t^ 
+ tu + v^) + s^(t + uf + stu{t^ + 3iii + V?) + t^u^{t + u)] , 

L = 2QV(s^ - 2t^) - 2QH{s^ - 2t^)[s{t - u) + t(t + u)] + Q'^s[s^(t'^ + u^) 
+ 2sH^{t + u) + st^{t^ + 6tu + u^) + At^u{t + u)], 

A^O. (A.6) 

^{T^) +9 ^ i^(i^) + cc pP(S^)] + g: 



277rM3(g2 + M|)2s'i(s + tfis + u)\t + uf' 

T = AQh'^ls^ - 2s\^ + s'^if - tu^) + s^{2i^ + M^u + 2^2-1^2 + ^^-^ + ^^(2^^ + g^s^ ^ ^^2^2 
^2tu^ + M^) + + + ^2)] _ 4g6^[s7(^2t - + sh{t + - 2s^(t^ + t^u + 2^^^ 

-M^) - sH{f + 6t2« + Itv" + 2m3) + - ?,iSi - leu" - 3t\^ + 2tu'' - u') + sH{t^ 
+At\i + 3t^u^ + 2t'^u^ + 3t?/ + u^) + + n)2(2t^ + 5tu + m^) + 2t^u\t + 

+2Q^t[2s*^(t - 2u) + 2/(t2 - M^) _ 56(3^3 _ 2t^u - tu^ - Su^) - 2s\2t^ + At^u + t^u^ 
-itu^ - 3u^) + s^{t^ - 6t% - 6t^u^ + 2fu^ + ?,tu^ - 4?/) + 2s^{t + M)^(2t^ + 2t^u 
-t^u^ + t?/ - 2m^) + shit + M)2(2t^ + 12t^?7, + 9t^M^ - 4tM^ - 2u^) + 2st\(t + uf(2t^ 
+Atu - u^) + 2t^u^{t + + 2g2stM[2s® + 2s\t + u)- s^{?,t^ + 2tu + S?;^) - 28^(2^3 
+bt^u + etM^ + 3m-'^) + s^(t^ - 4t^M - U^u^ - Atu^ + M^) + 2s-'^(2f' + If^u + llt^u^ 
+ m^ii^ + 6t?/ + 2?/) + s2(t + M)2(2t^ + I2t\ + ISt^-u^ + Qtu^ + 2m^) + 2stM(t + uf{2t^ 
+Atu + + 2t^u'^{t + u)% 

L = -2Qh^[s^ - 2s\^ + 2s^{t^ - tu^) + s2(4t^ + U^u + 4t2M2 + m^) + 2st(2t^ + 6t^M + bt^u^ 
+2tv^ + ?/) + 2t^{t + + u^)] + 2g^t[s^(3t - + 3sh{t + u) ~ 2s'\2t^ + 3t^u 

+Atu^ - u^) - 2s^t(t3 + St^M + Qtu^ + 2m3) + s^{<6t^ - 2t^u - m^u^ - 8t^u^ + 5tu^ - u^) 
+sH{Af + 12t\ + m^u^ + At\^ + 7tM^ + u^) + 2st\(t + uf{2t^ + 5tM + m^) + At^u^{t 
+uf] - Q^[s^{7t^ - Atu + u^) + 2sh{7t^ + Stw - 2^^) - s^(5t^ - 6t\ - t^u^ - I2tu^ 
+2u^) - 8sH{3t -u){t + uf - s^(8t^ + 76t^u + llOt^w^ + 40t^M^ - St^w^ + Stw^ - m^) 
+2s^t(4t^ + 2t^M - 22t\^ - 22t^u^ - fu^ ~~ ttf - 2u^) + sH'^{t + uf{At^ + 24t\ + St^^ 
-4tii^ - u^) + 4st^'u(t + M)3(2t2 + Atu - u^) + At'^u^{t + m)^] + 2g2s(s + u)[s\2t^ - tu 
W) + 2sH'^{3t + -u) + s^At"^ + m^u + 6t^w^ + 3iw=^ - 2u'^) - sHiAt"^ - llfu'^ - 12tu^ 
-u^) - s\6t^ + ISt^u + 12tV - + 2fu'^ + 2tu^ - u^) - sH{t + u)\2t^ + 12t\ 
+6t^u^ + 2tu^ + u'^) - 2st^u{t + uf{2t + 3u) - 2t^u^(t + uf] - s'^{s + uf{s + t + 
{t^ + u^) + 2s^{t^ - u^) + s^it"^ + 2^3x4 + 2i2ii2 + u"^) + 2st^M(t + m) + t^u^it + 'u)^], 

A = 0. (A.7) 
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^{^) +9 ^ ^(^) + cc ^Pi^^ + g: 



277rM3(g2 + Mlfs^s + f)^{s + uY{t + uf ' 
T = -8gi°t^(s + + m)[3s-'^ + s2(7t + 10m) + s(5t^ + IStw + Om^) + + St^w + 7tu^ 
+2u^] - AQ^t^[s^ + 7s\t + u) + s^{2At^ + 48t?i + 22^^) + 2s^(17t'^ + &2t^u + mtu^ 
+15m^) + /(7t^ + 112t3M + 219t\2 + imtu^ + 17m^) - s^(30t^ + 37t^?i - 91f\^ 
-ISlt^M^ - 56tM^ - 3m^) - sH{t + M)2(31t^ + Q2t^u - IStu^ - Au^) - st{llt^ + 69t^M 
+144t^M^ + 123t\^ + 42t\^ + Qtv!' + m^) - + + + 17tM^ + lu^)] 

+2Q^t[s^(4t - 2m) + 2s^{m'^ + 9iM - 77/2) ^ ^^7(^^25t3 + 206t2'u + 2btu^ - AAu^) 
+s^(239t^ + mbt^u + 507t2u2 - Qiw^ - 60m^) + s^(221t^ + 1022t^M + 1410tV 
+602t2ii=* - 35iii^ - 34w^) + s^(9U^ + 689*^-^ + 1548tV + 1318tV + SGU^ii^ 
-13tii^ - Qu^) + + M)2(16i^ + + 357^2^2 + mu^ + 2m^) + ^^^^(lOi^ 

+2At^u + 57^^-^^ + 127t^'u^ + m'^u^ - tu^ - 16u^) + 2st^(t + uf{At^ + 2t% 
-17t\^ - 19fu^ - mu" - Au'') + 2t\t + uYit" - Atu'' - u^)] - 2Q^[si°(3i' 
-2tu + -u^) + s^{25t^ + 7t\ - lltu^ + 7u^) + s^(100t^ + l?>2t\ - 2hfu'^ - 22tu^ 
+23ii^) + s^(208i^ + bm^u + 2m^u^ - m'^u^ - btu" + Mu^) + s^(217i^ 
+935*^1^ + 1162^^14^ + Zm\^ - mt^u^ + 29tvP + 22,u^) + s^(87i^ + 688^^-^ 
+1545iV + 1320iV + Am^u^ + 59^^^^ + mtu^ + 7w^) - s^{t + w)2(22i^ - 119i^'u 
-396^^-^^ - im^v? - 72t\^ - ntu^ - u^) - sH{28t^ + 152fu + 2h7fu^ + imt^u^ 
-163^^-^^ - 2Alt^u^ - lbd>t^u^ - bltu' - Au^) - sH^(t + u)^(6t'^ + 52t^u + VNlt^u^ 
+125iV + Alt^u" - 5tu^ - 7u^) - 2st^u{t + u)^{3t^ + I2t^u + 8*^^ + 2u^) - 2t\^{t 
+uf] + 2Q^s[s^°{f + v^) + s^{8t^ + + 6tu'^ + 8u^) + s^(32t^ + 56t^u + A2t^u^ 
+h2tu^ + 30w^) + s\7bt' + 232t^M + 2m^u^ + 22lt^u^ + 174^1/,^ + 57?/) + 5^(105^*^ 
+hm^u + 888t^M2 ^ smt^u^ + 576^2^^ + 280tM^ + b7u^) + 2s^(43t^ + 308^^^ + 7m^u^ 
+1015t^M^ + 77At^u^ + 389^2^^ + 119tM^ + Ihu') + /(i + uf{2,8t^ + 332t^M + 799^^-^^ 
^-7bAt\^ + ^.m'^u* + mu^ + 8u^) + s^{7f + 134t^M + 7Ut\? + 1847t^M2 + 2546^^-^^ 
+2066t^?/ + 993f'^M^ + 2Q2tW + 30tM^ + u^) + 2sHu{t + uf{8f + m^u + 167tV 
+193t^'u^ + lllt^u^ + 29tu^ + 2m^) + stW{t + 'u)^(lli^ + 34^2^ + 28^^^ + 7v?) 

+2tV(t + M)% 

L = 8Q^°t^(s + t)(t + m)[3s3 + s2(7t + 10m) + s{ht^ + IStw + 9^^) +e + 5t^u + 7tM2 + 2u^] 
+2Qh^[s^ - 2/(t + m) - s^(21t2 + 42tM + 23^^) - 2s^(29t^ + 62t'^u + 55tu^ + 22u^) 
-s\98t^ + 218t^u + net^u'^ + 90tM^ + 33m^) - 2s%5At^ + im^u + im^^i"^ + 
+ 13t?/ + 5?/) - s'^it + ufiJQf^ + im^u + Athi'^ + 6tu^ + u^) - 2st{llf + 69t^u 
+ lAAt\'^ + 123t\^ + 42^2?/ + 6t?/ + ?/) - 2f'^(t + uf{f + 9t'^u + ITtu'^ + 7m-'^)] 
-2g'^t[s^(3t - -u) + s^{f + 3tM + 2^2) - 5^(35^2 + 47^^^ - 3tu'^ - 23u^) - s%lQlt^ 
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-33ii^) - s^{m^ + ^m^u + 648iV + 620iV + 313^^-^^ + ^btu^ - IQu^) - s^{t + uf 
(30t^ + 126i^'u + lUfu^ + U2fu^ + Ihtu" - u^) + sH{<ot^ - 22t% - ISOt^u^ - 228iV 
-2Alt\'^ - 16UV - A7tu^ - u^) + 2st^{t + uf{At^ + 2t% - Ut^"^ - - I2tu^ 

-4«5) + 2t\t + u)\e - Atu^ - u^)] + g^[s^°(5t' - 8tu - i?) + 2s\4t^ - Q^u - IStw' 



-bu-") - s''{32t^ + 48t^u + 59t'u^ + 54tM^ + 35m^) - 4s''(30t^ + 63t% + 45iV + lOt'^u^ 
+ntu'^ + 13u^) - s^{179f + 596t\ + 69At\^ + 228t'^/ - 120^^?/ - 24:tu^ + 35u^) 
-2s^(78t^ + 404t^'u + 756t^M^ + 600t\^ + lAt^u* - 132t^M^ - ?>ltu^ + bu') - s^{t + uf 
(86^6 + 516t\ + 768t%2 ^ S60t\^ - U2t\* - 32tu^ + u^) - AsH{7f + 89fu + 338fu^ 
+596t^u^ + 535t^M^ + 212t^M^ - t^?/ - 19tu^ - u^) - sh\t + u^iif + 92t^u + 338t^u^ 
+444t^M''^ + 245t\^ + 30tu^ - 13u^) - 4:St^u{t + u)\3t^ + 12t\ + 8tu^ + 2u^) - 4t%\t 



+uf] + 2Qh{s + t + u)%s^u{5t 



+57t'^u^ 



u) + s^(2f^ + m\ + 21tu^ + 7u^) + s^^f^ + 2At\ 
mu^ + 16?/) + s^(3t^ + 13t% + SSt-'^M^ + 79t\^ + SOtw^ + 16m^) - s'^{f + lli^-u 
^22t%^ + 4:t^u^ - m'^u^ - 22tu^ - 7u^) - s\f + 15t% + 70t^u^ + 128t%^ + llSt^w^ 
+37t^M^ - 5tu^ - u^) - shu(t + uf{At* + 27t\ + A9tW + 2Qtu^ - u^) - st^u^(t + uf 
(5t^ + I9t^u + 18t?/ + 6^3) - 2tV(t + m)^] - .s^{s + u){s + t + uf\s''{t^ + 4tM + u^) 
+s^{Af^ + igtSv, + 18tM^ + hu^) + 2s^{3t* + m^u + 33t\'^ + 2U?/ + 5?/) + 2s^(2t^ 



3/J-6 



22t^M 



98t^u^ 



156f\-^ 



+20t'^M + 56f^?/ + 58t\^ + 26tM^ + 5m^) + s\t 
+A2tu^ + 5m^) + s'^u{t + M)2(5t^ + 3Qt\ + 33^^^^ + IGtw^ + u^) + stM2(t + uf{bt^ 
+7tu + Au^) + tV(t + -u)^], 
0. 



116tV 



(A.8) 



+ g ^ + cc 



3p(l) 



+ 5: 



1357rM3(g2 + Mlfs^s + + uY{t + m)^ ' 
-24g^°t^(s + t){t + u)[s^ + s2(3t + 4m) + s(3/ 
+4g®t2[s^ + bs\t + u) + 2s^{8t^ + mu + 7u^) 
+s^(61t^ + 168t^u + ISlt'^u^ + 102tu^ + 29w^) - 



7tu + 5u^)+t^ + 3t^u + 3tu^ + 2u^] 
s'\28t^ + 84^2^ + 82tu^ + 26u^) 



s%112t^ + 365i^ii + 443tV + 241t 
+68tu^ + 17u^) + s'^{t + uf{103t'^ + 234:t^u + U6fu^ + 20tu^ + 4w^) + st{A3t^ + 229t' 
+498iV + 551iV + 298t\'^ + 66tu^ + 7u^) +f{t + uf{7t'^ + 29^^-^ + A7fu'^ + A7tu^ 
+4ii^)] - 2Q^t[s^(4t - 2u) + 2s^(8i^ + 3tu - hu^) + s\m^ + 9At^u - 7tu'^ - 28u^) 
+s^(151i^ + AA3t\ + 3l9tW - 2htu^ - h2u^) + s^(225i^ + 8A2t'^u + 1102iV + 518i^w' 
-mu^ - b8u^) + s^(247i^ + 1025*^-^ + 1648iV + I2m^u^ + 385^^^^ + lltu^ - 3Au^) 
+s^{t + uf{180t^ + 628t\ + 545^^1*^ + lOGfu^ + 38tu'^ - 8u^) + sH{50t'^ + 5m% 
+1725iV + 237UV + 1543iV + 483^^1*^ + 92tu^ + 8u^) - 2st^{t + 1^)^(6^^ - bU'^u 
-229t^u^ - 2bbt^u^ - mu^ - lQ>u^) - 2t^{t + uf{3t^ - 3t^u - 32t\^ - A9tu^ - llu"^)] 



u 
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+62tu^ + llu"^) - 5^(64^5 - 27t'^u - 58UV - SOSt'^u^ - SBBtu"^ - A2u^) - s%22lf 
+557t^u - 54tV - 1270tV - UlSt'^u'^ - A37tu^ - 33u^) - s^Sllf + llSOfu 
+1263tV - 400tV - 1694tV - 1215iV - 336tu^ - 15u^) - s^{t + uf(266t^ + 7A5t^u 
+3m^u'^ - A6At\^ - A08t\^ - Ultu^ - 3u^) - sH{128t^ + SAOfti + lOSSt^w^ + 1756^^1*^ 
+69t\'^ - 969t^u^ - 694t\6 - 219tu^ - 28u^) - sH'^{t + uf{26t^ + 2A0fSi + 5A0t\'^ 
+321^2^^ - lllt\^ - 81tu^ - 29u^) - 2sfu{t + 'u)^(17t^ + 73t\ + 76t\^ + lOtu^ + 2u^) 
-14tV(i + uf] + 2Q'^s[s^^{3t'^ + 8tu + 3m^) + 6s\3f + 10t\ + Wtu'^ + 3u^) + s^(54i^ 
+228t^M + 3A6fu'^ + 2mu^ + A8u^) + s^(93t^ + 526t\ + llOSt^w^ + 106lfu^ + AUUi^ 
+75m^) + s^(99t^ + 764t^M + 2132t^M2 + 2872t^M''^ + 1980t^M^ + Qh2tu^ + 7bu^) + s^(72t^ 
+734t^M + 2618t^M^ + 4630t^M^ + 4480t^M^ + 2370t^M^ + &l2tu^ + 48m^) + s^{t + uf{3&f^ 
+404t^M + 1297t^?i^ + 1774t^M^ + llA3thi'^ + 332tM^ + I8u^) + s^{9f + 188t^M + llA2t^u^ 
+325U^M^ + 5142t^M^ + 4834t^M5 + 2729t^M^ + 874t\^ + 128tM^ + 3u^) + 2sHu{t + uf{m^ 
+ 138t^i + 387tSi^ + A9hf^ii^ + 311t^?/ + imtu^ + IOm^) + sfu'^{t + uf{37t'^ + IQlt^u 
+23At^u^ + 127tM-^ + 29m'*) + Ui^u^it + u)% 
L = 24Q*°t^(s + t){t + + s^{3t + 4m) + s{3t^ + 7tu + Sm^) + + 3t^u + St?/^ + 2u^] 

~2Q^t^[s^ + U8\t + u) + s%55t^ + llOtu + 53u^) + 2s\67t^ + 162thj, + UUa? + 46m^) 
+s^(254t^ + 666f^n + 652tW + 322tM^ + 83m^) + 2s^(160t^ + 521t\ + 608t^u^ + 319t^u^ 
+91tu^ + 19u^) + s\t + M)2(230t^ + 516t^u + 304t\2 + 45^^3 ^ 7^4^^ ^ 2st(43t^ + 229^^^ 
+mt\^ + 551f^M''^ + 298t^M^ + 66tu^ + 7u^) + 2t^{t + M)^(7t^ + 29t^M + A7t\'^ + mu^ 
+4u^)] + 2Q%s^{3t -u) + s\51t^ + 37tu - Uu^) + s^(221t^ + 389t\ + mtu^ - 53u^) 
+s^(467t^ + 1269t^u + 1049t^M^ + 155tu^ - 92u^) + s^(610t^ + 2090t^M + 2674:t^u^ 
+ 1434t\3 + 161tu^ - 83u^) + /(542t^ + 2196t^u + 3484^^^^ ^ 2782^^^ + 1115t\^ 
+125tu^ - 38u^) + s%t + uf{302t^ + 990t% + 870t\^ + 3A2t\^ + 67tu'^ - 7u^) + 5^^(70^^ 
+Q>7Qt^u + 2020^^-^2 + 2756tV + 1859tV + U3t\^ + \27tu^ + 7ii^) - 2st^{t + 'u)2(6i^ 
-Set^M - 229iV - 255^^1*^ - mu^ - l&vP) - 2f{t + M)^(3t^ - 3fu - 32t\'^ - mu^ 
-llu^)] - Q^[s^%13t'^ + 8tu + 7u'^) + 2s'^(74i^ + m'^u + tu^ + 19u^) + s^{b88t^ + 88Qt^u 
+61i^ii^ - imu^ + 89u^) + 4s^(298i^ + 7hQt% + 413^^1*^ - IGli^ii^ - 93tu'' + 29ii^) 
+s^(1421i^ + h2m^u + 5998iV + 1312i^ii^ - 1504i^w^ - 536iw^ + 89u^) + 2s^(536t^ 
+2616t^'u + 4606i^w^ + 3264iV + 232iV - 73Afu^ - 2mtu^ + I9u') + + uf{b3Af 
+2328t^u + 2692^^1^^ + 89m^u^ - 37m'^u^ - 18Atu^ + 7u^) + 4s^i(43i^ + 397^^-^ + imt^u^ 
+1596t^u^ + lOeOtV + 279iV - 55^^-^^ - 5Atu^ - 7u^) + sH'^{t + uf{28t^ + 444i^'u 
+1178iV + 920iV + 103i^M^ - 2mu^ - A3u^) + Ast^u{t + 'u)^(17i^ + 73^^-1* + 7m'^u'^ 
+mu^ + 2ii^) + 28tV(t + -u)^] + 2Q'^s[s^^{8f + llii* + 7ii^) + s^{m^ + 78t^'u + Qhtu^ 
+Ahu^) + s^(177i^ + 287t\ + 162^^x4^ + 187^^^ + 127ii^) + s^(295i^ + Qhbt% + 232iV 
-2Qt^u^ + 3l3tu^ + 205w^) + s^(290i^ + 92bt^u + 438i^w^ - 844^^-1^3 _ 540^2^^^ ^ ggg^^s 
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+53tS - m^u^ - - mit^u^ - ^m'^u^ - ntu^ + lu') - sV{t + ufi^bt"^ + 2iAt^u 

+227t\'^ + 99tu^ + u^) - 3sfu^(t + uf{lt^ + m'^u + 22tu^ + Q>u^) - 14iV(i + uf] 
-s^{s + ii)(s + t + uf[s\7t^ + mu + 7u^) + s^(28t^ + Aht^u + mu^ + 31m^) + 2s^(2W^ 
+33t^M + 51t\2 + mtu^ + 29m^) + 2s^(14t^ + 24t^M + 34^^^ + m'^u^ + 94tM^ + 29m^) 
+s2(7t^ + 18f'^?i - 2t^M2 + mt\^ + ISOt^M^ + l2>Stu^ + 31?i^) + s^'u(i + uf{2,t'^ - lAt\ 
+23^2-^2 + mtu^ + 7ii^) + ?,stu'^{t + + 3iw + Au^) + 7i2'u^(i + u)% 
0. (A.9) 



10(8) 

"^0 



9- 



4 2 



127rM(Q2 + M^ysH{s + t)2(s + M)2(t + m)2 ' 
T = 9t;^g^(g2t + SM)[2g^t2M2 + 2Q^stu{st + su + t^ + tu + u^) + 

+ 2s^(t + m) + 3s^(t + + 2s^{t + m)^ + + 2t^u + 3^^^ + 2tv? + -u^)] 
+ 10a^Q^s^t2(Q2t + su){s + t + uf, 
L = -9vIqH[2QH'^u^ + 2Q'^stu{st + su + t'^ + tu + 2m^) + /(t + m)^ 
+ 2s^{t^ + 2tSi + 2tu'^ + + s2(t^ + 2t^u + 3t^M^ + 2tu^ + 2m^)] 
- balsH{QH + sm)2(s + t + uf, 
A = Q. (A.IO) 

1/(17) + ^ ^ v{v) + cc [^^f ^] + g: 

4 2 



367rM(Q2 + MlfsH{s + + + m)2 ' 
T = 10t;fg2t{2g6t2(g2 ^ ^2^) _ q4^[s3(^3^ _ + 352^^^ + + 2sf{t - u) + 2t^{t + u)] 

+ Q^s[s^{t - uf - 2s'^tu{t + u)- sfu{2t - n) - 2f^;i(t + u)] - s^{f + tu + u^) 

- s^{t + uf - sHu{f + 3tu + u'^) - sH\'^{t + u)} + 9alQ^{2QH^u'^ - 2QH'^u[s^ 

- s^{t + u) + s{t^ + tu- 2u'^) + 2tu{t + u)] - Q^st[s^ + 2s^(t + u) + s^{f + 2tu 
+ 3m^) + 2sH{t'^ + 3tu + 2u^) + s(t^ + 6fu + Ibt^u'^ + 12tu^ - u"^) + Atu'^{t'^ + 3tu 
+ 2u^)] - s'^u - 2s%{t + u)- s^{2t^ + bfu + Uu^ + Zu^) - 2s\t'^ + At\ + ifu^ 
+ htu^ + u^) - s^iZt"^ + I2t\ + nfu^ + mu^ + u^) - 2sV{t^ + Zfu + Atu^ 

L = 5vIQH{2QH''{s'' - 2^) - 2QH{s'' - 2e)[s{t - u) + t{t + u)] + s\e + u"") 

+ 2sV(i + li) + sH'^it^ + mu + u^) + Ast^u{t + u)}- 9al{2Qh^u'^ + 2QH'^u[s^ 
+ s^{t + u)- s{t^ + tu- 2u^) - 2tu{t + u)] + Q'^stis^ + 2s^{t + u) + s^u{t + 5u) 

- s^{2t^ + 7t\ + tu^ - Au^) - sit"^ + 6t^u + lAt^u^ + %tu^ - 3w^) - Atu'^{t^ + ?,tu 
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+ 2u'')] - Qh^[s\t -u) + s\3t'' + tu- 2u^) + s^{3t^ + St^u + Uu^ - 3u^) + s^t^ 
+ 9t^u + m'^u^ + 9>tu^ - 2u^) + suiZt"^ + m^u + 2Qt^u^ + 9tu^ - u^) + tu^{t 
+ uf{t + bu)] - s^u{s + t + + s{t + li) + u{t + li)]^}, 
A = 0. (A.ll) 



T = 9v^Q2{8g^t^[2s^t - 2s\'^ + sH^{2t - Zu) - sHu{t^ + 3tu - 3u'^) - stu{2t^ + t'^u 

- tu^ - V?) + fv?{2t + u)] - 2Qh[As\t + u) + As^{ht^ + 2u^) + As^{bt^ + 

- 2tM2 + 4^=^) + 2s^(9t^ - 9t^M + 4t2'u2 _ gt^s + 4m^) + s^(12t^ - 2t^'u - 29^3^2 
+ SSt^?/"^ - I2tu^ + 4m^) - s^t(2t' + 20t^?7, + 25t\^ - t^u^ - IQiu^ + Au^) - st\ 
X (12t^ + 2Qt^u + 2t\2 + ^ ^4^) _ ^3^^2(^2t3 _ Qfu - 15tM^ - 7u^)] + 2Q'^[2s^ 

X (3t^ + Ui - 2u^) + s\2lt^ + 3t\, + lOtM^ - ^u^) + s^(28t^ + Ibt^u - 15t\^ + 22^^^ 

- I2u^) + s^(27t^ + llt\ + 8t^M^ - 20t^M^ + 2Qtu^ - 8?/) + s^(16t^ + 17t^?i - 20t^M2 
+ 2t^u^ - 33t^?/ + 18t?/ - 4?/) + s3t(2t^ + 8t^M - 3^^^^ - m^u^ - m'^u'^ - ITtw^ 

+ 6m^) - s2t^(2t*^ + lAt'u + m\i'^ + m^a^ + A&t^u^ + 2ltu^ + Au^) - 2st^u{2t> 
+ lOt^M + m^^i"^ + 20t^?7;'^ + 19t?/ + 7m^) - 2t^M^(t + uf] - Q'^s[s^{yt^ - 5tu - 12m^) 
+ s^(25t^ + 3t^u - IStu'^ - 36u^) + s^(37t^ + 25f'^M - Ut'^u^ - 20tu^ - 60m^) + s^(39t^ 
+ 39t^?7, + mt^u^ - 18t\^ - 6tu'^ - 60?/) + s^(29t^ + 83t^M + 72t'^u'^ + 88t\^ + AOt^^ 
+ 22t?7;'' - 36?/) + s^(9t^ + 75t^?7, + 148/?/ + 176/?/ + 178/?/ + 102/?/,^ + 22t?/ 

- 12u^) + sHu{22t^ + 107/?i + 199/??^ + 211/?^^ + 177/??^ + 73t?i^ + 9??^) + sfu^{17t^ 
+ 69/?/ + 107/?/^ + 105/?/^ + 71tu'^ + 21u^) + At^u^{t + uf] - s^{s + u)[7s\{t + ??) 

+ /??(25t2 + 38t?? + 21??2) + 5^(2/ + 47/?? + 88t\'^ + 78tu^ + 35??^) + s^(4t^ + 63/?? 
+ 132/?/^ + 156t\;'^ + 98t?i^ + 3bv^) + s^{2f^ + 47/?i + 136/?/ + 190/?/ + 156t2??^ 
+ 78t??^ + 21??^) + 5^^(13/ + 70/?? + 136/u^ + 132/??^ + 88/?/ + 38tu^ + 7u^) + stu^ 
X (13/ + 47/-^ + m^u^ + 47/?x=^ + 2btu^ + lu^) + 2t\^{t + uf]} + IQalQHiAQH"^ 
X [2s^ + 5^(5* + 7??) + s(4/ + lltu + lu^) + {t + uf{t + 2??)] + 2QH\2s^ + s^{7t 
+ 9u) + 5^(5/ + 27tu + I2u^) - 5^(13/ - 7/w - ZAtu^ - Au^) - 5^(23/ + h7t\ 
+ - 14iii^ + 2??^) - s(12/ + 52/-^ + 73/^2 + 33fu^ + tu'^ + -u^) - i(2/ 

+ 12/-^ + 27t\^ + 29/-^^ + mu'^ + w^)] - QH[2s\bt - 2u) + s^(4l/ + 37tu - ISw^) 
+ 5^(56/ + 147/-^ + Ahtu^ - 2Au^) + s^(19/ + 164/ii + 194/w^ + mu^ - 8w^) 

- 5=^(16/ - 12/-^ - 154/1^2 - 117/w=^ + btu"^ - Au^) - s'^{8t^ + 66/-^ + 95/1^^ 
+ 8/w^ - 17fu^ + lOtu^ - 2u^) + 2st(2/ - 8/w - 53/ii^ - 82/ii^ - A5fu^ 

- 8tu^ - 2u^) + 2/(i + 'u)2(/ - 8/-^^ - 12tu^ - 3u'^)] + Q'^[s^(t + u) + s^(ll/ 
+ tu + 6u'^) + 5^(40/ + 36t\ - 3tu'^ + 17u^) + s^(56/ + 144/-^ + 57/ii^ + Qtu^ 
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- 182iV - 228iV - 117^^1^^ - b^tu^ - Q>u^) - s^(30i^ + 126iS + llbt^u^ 

- 81iV - 158iV - m'^u^ - mu^ - u^) - sH{8f + 72t% + 191i V + ISlt'^u^ 
+ 16tV - 60t\^ - 2mu^ - 7u^) - 2sfu{t + uf{ht'^ + 2lt^u + 2lt\'^ + 2tu^ 

- u^) - At^u^{t + uf] - s{s + u)[s^{t + u) + Qs\t^ + tu + u^) + /(19t^ + 2'dt^u 
+ 27tu^ + I7u^) + 2s^(17t^ + A7t^u + SOt^?;,^ + SSt?/'^ + I2u^) + /(35t^ + im% 
+ 254t3M2 + 202^2-^3 + 92i^4 ^ -^^^5^, ^ 53^^20^^ + 145^^1* + 2>m\^ + 370t^M^ 

+ 216^21^^ + Qhtu^ + 6m^) + s2(5t^ + m^u + 217i^ii^ + 346t^M^ + 286^^^^ + 125t\^ 
+ 2Qtu^ + ?/) + stu{t + M)2(9t^ + 40t3M + bQt^u^ + 28tM3 + hu^) + 4t2M2(t + li)^]}, 
lSvlQ*{AQ^t^[sH - - 2s-'^f'^ + sh'^u{t + 3m) + stu{2t^ + t\ - tw^ - u^) 

- t\\2t + u)] - QH[2s'^{t + u)+ 4/(2t2 + u^) - s^{t^ - t\ + 2tu^ - 8u^) - /(21t^ 
+ 13t\ + 18fu^ - 2tu^ - Au^) - 2s^{<6t^ + lOtS + t^u^ + At^v!^ - 6t?/ - + 2s2t 

X {t^ + 10t% + llt^M^ ^ 9^2^3 ^ g^^4 ^ 5^5) ^ st2M(12t^ + 26^^^ + 2t^U^ + tM^ + 

+ t3«2(2t3 _ 6^2^ _ ^5^^2 _ 7^^3)] ^ 2Q^[s\e - u") + /(3t3 + tu" - 2^3) _ s\2t^ 
+ 2t^M + 2t\2 - tu^ + 3m^) - f^\l2t' + 16t^M + 16t^M^ - t^u^ + to^ + 2u^) - /(lOi^ 
+ 28t^M + 28tS,2 + 13t^M-'^ - 7t^u^ + 3t?/ + m^) - sH{f + 9t^M + IStS^ + t^w^ 

- ISt^M"^ - mAi^ + 2?/) + s^f^^it"^ + 7f^M + 13t^?i^ + m^u^ + 35t2M^ + 2Atu^ + 6m^) 
+ sf\,(2t^ + m\i + 16f^7,^ + 2Qt^u^ + 19tM^ + 7m^) + f'u^it + m)^] + s{s + M)[2s^M(t 
+ m) + 2s^M(3t2 + 3tM + 2^2) + 5^(5^4 + 15^3^ + i8tV + lAtu^ + 6m^) + s^{lht^ 

+ 38t^M + 53t^M2 + 40t\^ + 22tM^ + Au^) + s^(15t^ + 52t^M + 88t^M^ + m^u^ + Alt^u^ 
+ mu^ + 2m^) + sh{U^ + 32t^M + 7U\^ + 90t^M-'^ + 68t^M^ + 34tM^ + 9m^) + st^uiyt' 
+ 31t^M + 47tV + 39tV + 23tM^ + 7u^) + 2t^u^{t + w,)^]} - lQalt{AQ^h^[2s^ 
+ s^(5t + 7m) + s(4t2 + UtM + 7u^) + (t + M)2(t + 2m)] - 2g*^t2[2s'^ + 2s^(4t + 5m) 
+ 2s^(9t2 + lltu + 9^2) + s^(29t^ + A2>t^u + 15^7/2 + 15-^3) + ^2(^27^^ + 73t3M + Alt^u^ 
+ tM^ + 6m^) + s(12t^ + 52t^M + 72,t^u^ + 33^^^^ + tu^ + m^) + t(2t^ + 12t^M + 27tV 
+ 29i^M^ + 13iM^ + M^)] + QH[s\m - Au) + 5s^(9i^ + 8tM - 4^^) + s\m^ + 153t^M 
+ 73^^^ - 36^3) + s^(87i^ + 22lt^u + 226i2M2 + mtu^ - ZQu^) + s'^(56i^ + im'^u 
+ 231iV + l%Zfu^ + 82iM^ - 12m^) + s^{lAt^ + 106t^M + 201i^M^ + 145i^M^ + 69i^M^ 
+ mvP - 2u^) - 2st{2t^ - St^u - 53iV - 82t^u^ - Abfu^ - Stu^ - 2u^) - 2f{t + uf 
X (i^ - ^fu'' - mu^ - 3m^)] - Q^[s\t + m) + s\lAe -tu + 6m') + s''{A8t^ + 35t\ 

- 18tu^ + 16m^) + s^(78i^ + 148t^M + 17i'M' - 54tM^ + 22m^) + s\75t^ + 2A0t% + 214tV 

- 30i'M^ - 86iM^ + 16m^) + s\A8t^ + 229t^u + 360t%'^ + 216t\^ - ASt'^u'^ - 7Atu^ 

+ 6m^) + s^{20f + 152t% + 339t^u^ + 326iV + 134iV - 17i'M^ - 31iM^ + m^) + sH{Af 
+ 62t% + 210t^M' + 275iV + 150tV + 28i'M^ - Atu^ - 5m^) + 2st\{t + uf{5t^ + 21t^u 
+ 2ltW + 2tu^ - u^) + AtW{t + m)^] + Q'^s[2s^{t + m) + s^{l2f + 12iM + 13m') + s\28t^ 
+ 41i'M + 37iM' + 37m=^) + s^{m^ + 73i=^M + 50i'M' + 62tM^ + 60m^) + 5^(18^^ + 59i^M 
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+ 37ii^) - s^u{m^ + 84^5-^ + 225iV + 279^^-^=^ + 162^^1^^ + 21tu^ - 13ii^) - s^u{t + uf 
X {At^ + AQt% + m^u^ + hlfu^ + mu^ - 2u^) - 2stv?{t + uf{At^ + Ihfu + 9*^^ + u^) 
- At^u^{t + uf] - s^{s + + i + + -u) + s^(3t^ + 4iii + Au^) + s^(3t=* + 7t^it 

+ lltu^ + 7ii^) + + 6i^ii + 13i^ii^ + Utu^ + 7ii^) + s^u{2t'^ + St^ii + 13t^'u^ + lltu^ 
+ 4x4^) + sv?{t + 1^)^(2^^ + 2tii + u^) + + uf]}, 

0. (A.12) 
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Figure 1: Representative Feynman diagrams for the partonic subprocesses u + a u + 
cc[n] + a, where a = q,q,g and n = '^Si^\'^Pj^\ '^SQ^\ '^S[^\'^Pf'\ There are six diagrams 
of the type shown in part (a), two ones of the type shown in part (b), and two ones of 
the type shown in part (c). There arc two more diagrams that are obtained from the 
diagrams of the type shown in part (b) by replacing the external gluon lines with quark 
ones. The CS process only proceeds through the diagrams shown in part (a). 
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Figure 3: LO NRQCD (upper histogram) and CSM (lower histogram) predictions for 
the distribution of prompt J/ip inclusive production in nondiffractive vN NC DIS 
appropriate for the CHORUS and NOMAD experiments. 



32 



-2 
10 




y 

Figure 4: LO NRQCD (upper histogram) and CSM (lower histogram) predictions for the y 
distribution of prompt J/ip inclusive production in nondiffractive i/N NC DIS appropriate 
for the CHORUS and NOMAD experiments. 
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Figure 5: The total cross section oiu + N^u + J/ilj + X measured by CHORUS [|T2 
(EX) is compared with the predicted diffractive (DI), CS, and CO contributions. 
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